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Abstract 

Strong electron-phonon interactions in cuprates and other high-temperature superconductors 
have gathered support over the last decade in a large number of experiments. Here I briefly in- 
troduce the Frohlich-Coulomb multi-polaron model of high-temperature superconductivity, which 
includes strong on-site repulsive correlations and the long-range Coulomb and electron-phonon 
(e-ph) interactions. Our extension of the BCS theory to the strong-coupling regime with the 
long-range unscreened electron-phonon interaction could naturally explain isotope effects, uncon- 
ventional thermomagnetic transport, and checkerboard modulations of the tunnelling density of 
states in cuprates. 
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I. INTRODUCTION: THE FROHLICH-COULOMB MODEL 

A significant fraction of research in the field of high-temperature superconductivity 
suggests that the interaction in novel superconductors is essentially repulsive and unretarted, 
but it provides high without any phonons. A motivation for this concept can be found 
in the earlier work by Kohn and Luttinger 0], who showed that the Cooper pairing of 
repulsive fermions is possible. However, the same work showed that of repulsive fermions 
is extremely low, well below the mK scale. Nevertheless, the BCS and BCS-like theories 
(including the Kohn-Luttinger consideration) heavily rely on the Fermi-liquid model of the 
normal state. This model fails in cuprates, so that there are no obvious a priory reasons 
to discard the dogma, if the normal state is not the Fermi-liquid. Strong onsite repulsive 
correlations (Hubbard U) are essential in undoped (parent) cuprates, which are insulators 
with the optical gap about 2eV or so. Indeed, if repulsive correlations are weak, one would 
expect a metallic behaviour of a half-filled d-hand of copper, or, at most, a much smaller 
gap caused by lattice and spin distortions (i.e. charge and/or spin density waves [4]). It is a 
strong onsite repulsion of (i-electrons which results in the "Mott" insulating state of parent 
cuprates. Different from conventional band-structure insulators with completely filled or 



illed band as 
\5\. 



empty Bloch bands, the Mott insulator arises from a potentially metallic hal 
a result of the Coulomb blockade of electron tunnelling to neighboring sites 

However, independent of any experimental evidence the Hubbard U (or t — J) model 
shares an inherent difficulty in determining the order. While some groups claimed that it 
describes high- Tc superconductivity at finite doping, other authors could not find any super- 

n 

conducting instability without an additional (i.e. e-ph) interaction [6]. Therefore it has been 
concluded that models of this kind are highly confiicting and confuse the issue by exagger- 
ating the magnetism rather than clarifying it There is another serious problem with the 
Hubbard-U approach to high temperature superconductivity in cuprates. The characteristic 
(magnetic) interaction, which might be responsible for the pairing in the Hubbard model, is 
the spin exchange interaction, J = At'^/U, of the order of 0.1 eV. It turns out much smaller 
than the (inter-site) Coulomb repulsion and the unscreened long-range (Frohlich) e-ph inter- 
action each of the order of 1 eV, routinely neglected within the approach. There is virtually 
no screening of e-ph interactions with c— axis polarized optical phonons in cuprates because 
the upper limit for the out-of-plane plasmon frequency (< 200 cm-i)^ is well below the 
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characteristic phonon frequency, uo ^ 400 - 1000 cm ~^ . As a result of poor screening, the 
magnetic interaction remains small compared with the Frohlich interaction at any doping. 
Hence, any realistic approach to superconductivity and heterogeneity in cuprates should 



treat the Coulomb and unscreened e-ph interactions on an equa^ 



footing. 



We have developed a so-called "Frohlich- Coulomb" model 0, 1^ to deal with the strong 
Coulomb and long-range e-ph interactions in cuprates and other doped oxides. The model 
Hamiltonian explicitly includes the long-range electron-phonon and Coulomb interactions 
as well as the kinetic and deformation energies. The implicitly present large Hubbard term 
prohibits double occupancy and removes the need to distinguish fermionic spins. Introducing 
spinless fermionic, Cn, and phononic, d^a, operators the Hamiltonian of the model is written 
as 



- "^Uagaim - n)(e„ ■ Uin_n)c|^Cn(rf|^„ + c/^a) + "^1^^ HLa'^ma + ^ j , (1) 
nm ma ^ ' 

where Sq, is the polarization vector of ath vibration coordinate, Um-n = (m — n)/|m — n| is 
the unit vector in the direction from electron n to the ion m, (70(111 — n) is a dimensionless 
e-ph coupling function, and K(n — n') is the inter-site Coulomb repulsion. (70,(111 — n) is 
proportional to the force acting between an electron on site n and an ion m. For simplicity, 
we assume that all the phonon modes are non-dispersive with the frequency cUa- We also 
use h = 1. 

The Hamiltonian, Eq.(l), can be solved analytically in the extreme case of a strong e-ph 
interaction where the e-ph dimensionless coupling constant is large, A = Ep/zt > 1, by using 
1/A multipolaron expansion technique Here Ep = Xlno 'Un)^, is the polaron 
level shift about 1 eV and zt is the half-bandwidth in a rigid lattice. 

The model shows a reach phase diagram depending on the ratio of the inter-site Coulomb 
repulsion Vc and the polaron (Franc-Condon) level shift Ep pij]. The ground state is a 
polaronic Fermi liquid at large Coulomb repulsions, a bipolaronic high-temperature super- 
conductor at intermediate Coulomb repulsions, and a charge-segregated insulator at weak 
repulsion. The model predicts superlight bipolarons with a remarkably high superconduct- 
ing critical temperature. It describes many other properties of cuprates jof, in particular 
band-structure isotope effects, normal state thermomagnetic transport and real-space mod- 
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ulations of the single-particle density of states (DOS) as discussed below. 



II. BAND-STRUCTURE ISOTOPE EFFECT 

The isotope substitution, where an ion mass M is varied without any change of electronic 
configurations, is a powerful tool in testing the origin of electron correlations in solids. 
In particular, a finite value of the isotope exponent a = —d\nTc/d\nM proved that the 
superconducting phase transition at T = Tc is driven by the electron-phonon interaction 



in conventional low-temperature superconductors Advances in fabrication of isotope 

substituted samples made it possible to measure a sizable isotope effect also in many high- 
temperature superconductors. This led to a general conclusion that phonons are relevant 
for high Tc. Essential features of the isotope effect on Tc, in particular its large values in 
low Tc cuprates, an overall trend to decrease as Tc increases, and a small or even negative a 
in some high Tc cuprates, could be understood in the framework of the bipolaron theory of 
high-temperature superconductivity [l^ . 

The most compelling evidence for (bi)polaronic carries in novel superconductors was 
provided by the discovery of a substantial isotope effect on the (super) carrier mass [isl 



providec 

Q,li3- 



The effect was observed by measuring the magnetic field penetration depth Xh 
of isotope-substituted copper oxides. The carrier density is unchanged with the isotope 
substitution of O^^ by O^*^, so that the isotope effect on Xh measures directly the isotope 
effect on the carrier mass m*. A carrier mass isotope exponent = dlnm* /dlnM was 
observed, as predicted by the bipolaron theory In ordinary metals, where the Migdal 
adiabatic approximation is believed to be valid, am = is expected. However, when 
the e-ph interaction is sufficiently strong and electrons form polarons, their effective mass 
m* depends on M as m* = mexp{'yEp/uj). Here m is the band mass in the absence of 
the electron-phonon interaction, 7 is a numerical constant less than 1 that depends on the 
radius of the electron-phonon interaction, and is a characteristic phonon frequency. In 
the expression for m*, only the phonon frequency depends on the ion mass. Thus there is 
a large isotope effect on the carrier mass in (bi)polaronic conductors, am = (1/2) ln(m*/m) 
[1^ . in contrast with the zero isotope effect in ordinary metals. 



Recent high resolution angle resolved photoemission spectroscopy (ARPES) 17[ provided 
another compelling evidence for a strong e-ph interaction in cuprates. It revealed a fine 
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phonon structure in the electron self-energy of the underdoped La?_.rSrrCu04 samples jisl ]. 
and a complicated isotope effect on the electron spectral function in Bi2212 [l9|. Polaronic 
carriers were also observed in colossal magneto-resistance manganite including their low- 
temperature ferromagnet ic p hase, where the isotope effect on the residual resistivity was 
measured and explained 20|. 

These and many other experimental and theoretical observations point towards unusual 
e-ph interactions in complex oxides, which remain to be quantitatively addressed. We have 



performed quantum Monte Carlo (QMC) simulations j2l| of a single-polaron problem with 
the long-range e-ph interaction, Eq.(l), in the most relevant intermediate region of the 
coupling strength, A ~ 1, and of the adiabatic ratio, uo/t ~ 1, where any analytical or 
even semi-analytical approximation (i.e. dynamic mean-field approach in finite dimensions) 
may fail. More recently we have calculated the isotope effect on the whole polaron band 
dispersion, et, and DOS of a two-dimensional lattice polaron with short- and long-range 
e-ph interactions by applying the continuous-time QMC algorithm [2^. Unlike the strong- 
coupling limit [3] the isotope effect depends on the wave vector in the intermediate region 
of parameters. It also depends on the radius of the e-ph interaction. If we define a band- 
structure isotope exponent as 

it does not depend on the wave- vector k in the extreme strong-coupling limit, A ^ 1. E-ph 



interactions do not change the band topology in this limit [12|, so that ab is the same as 
am- Notwithstanding QMC results show that the isotope exponent becomes a nontrivial 
function of the wave vector, = ^^(k), in the intermediate-coupling regime, because e-ph 
interactions substantially modify the band topology in this regime. 

The isotope exponent ab(k) are presented in Figs. (1,2), for the small Holstein polaron 
(SHP) with the short-range interaction, and for the small Frohlich polaron (SEP) with the 



long-range interaction [2^. The polaron spectra are calculated for two phonon frequencies, 
uj = 0.70 1 and uj = 0.66 1, whose difference corresponds to a substitution of O^^ by O^^ in 
cuprates. There is a significant change in the dispersion law (topology) of SHP, Eig.l, which 
is less significant for SEP, Eig.2, rather than a simple band-narrowing. As a result, the iso- 
tope exponent, ab(k) ^ 8[e^ is k-dependent, Eig.l and Eig.2 (lower panels). The 
strongest dispersion of ab is observed for SHP. Importantly, the isotope effect is suppressed 
near the band edge in contrast with the k-independent strong-coupling isotope effect. It is 



less dispersive for SFP, especially at larger A. 

The coherent motion of small polarons leads to metallic conduction at low temperatures. 
Our results, Figs. 1,2 show that there should be anomalous isotope effects on low-frequency 
kinetics and thermodynamics of polaronic conductors which depend on the position of the 
Fermi level in the polaron band. Actually, such effects have been observed in ferromagnetic 
oxides at low temperatures [2^, and in cuprates [3,13, 15 1. To address ARPES isotope 
exponents one has to calculate the electron spectral function A(k, E) taking into account 
phonon side-bands (i.e. off diagonal transitions) along with the coherent polaron motion 
(diagonal transitions). Using the 1/A expansion and a single phonon mode approximation 



one obtains 



where 



00 

A{k, E) = J2 (k, E) + Ai"-^ (k, E)] , (3) 
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and 



A^"Hk,i^) = ^-g^+/"^ E i{\lM\%E-lu-i{^.i)]. (4) 

' qi.-qi'"=i 

Here Z = exp{—Ep/uj), C,(k) = — fi, fi is the chemical potential, n{E) = [exp(/3_E) + 

k^ = k ± X]t=i and 7(q) is the Fourier transform of the force, 7(q) = V^J2n9(^)(^ ' 

Un) exp(iq ■ n). 

Clearly, Eq. (jH)) is in the form of a perturbation multi-phonon expansion. Each contri- 
bution A;^'' (k, E) to the spectral function describes the transition from the initial state k 
of the polaron band to the final state k^'' with the emission (or absorption) of / phonons. 
Different from the conventional Migdal self-energy the electron spectral function com- 
prises two parts in the strong-coupling limit. The first (/ = 0) k-dependent term arises from 
the coherent polaron tunnelling, Acoh(k, E) = A[,~^(k, E) + A[,"*"^(k, E) = Z6{E - Ck) with 
a suppressed spectral weight Z 1. The second incoherent part AincohO^,E) comprises 
all terms with Z > 1. It describes excitations accompanied by emission and absorption of 
phonons. We note that its spectral density spreads over a wide energy range of about twice 
the polaron level shift Ep, which might be larger than the unrenormalized bandwidth 2zt 



in the rigid lattice. On the contrary, the coherent part shows a dispersion only in the en- 
ergy window of the order of the polaron bandwidth. It is important that the incoherent 
background AincohiM, E) could be dispersive (i.e. k-dependent) for the long-range inter- 
action. Only in the Holstein model with the short-range nondispersive e-ph interaction, 
7(q) = constant, the incoherent part has no dispersion. 

Using Eq. one readily predicts the isotope effect on the coherent part dispersion and 
its spectral weight Z, and also on the incoherent background because Z, 7(q), and uj all de- 
pend on M. While our prediction is qualitatively robust it is difficult to quantify the ARPES 
isotope effect in the intermediate region of parameters. The spectral function, Eq. (jH}, is 
applied in the strong-coupling limit A ^ 1. While the main sum rule J^^dEACk, E) = 1 
is satisfied, the higher- momentum integrals dEEPA{]<., E) with p > 0, calculated with 
Eq. differ from their exact values j23| by an amount proportional to 1/A. The differ- 
ence is due to partial "undressing" of high-energy excitations in phonon side-bands, which 
is beyond the lowest order 1/A expansion. The role of electronic correlations should be 
also addressed in connection with ARPES. While the results shown in Figs. 1,2 describe 
band-structure isotope effects in slightly-doped conventional and Mott-Hubbard insulators 
with a few carriers, their spectral properties could be significantly modified by the polaron- 
polaron interactions [2^, including the bipolaron formation 25] at finite doping. On the 
experimental side, separation of the coherent and incoherent parts in ARPES remains rather 
controversial. 



III. THERMOMAGNETIC TRANSPORT 

There is much evidence for the crossover regime at T* and normal state charge and 
spin gaps in the cuprates 0, |2^. Within the Frohlich- Coulomb model these energy gaps 
could be understood as being half of the binding energy, Ap, and the singlet-triplet gap 
of preformed bipolarons, respectively Notwithstanding some "direct" evidence for the 
existence of a charge 2e Bose liquid in the normal state of cuprates is highly desirable. 



Mott and Alexandrov 27| discussed the thermal conductivity k; the contribution from the 
carriers given by the Wiedemann-Franz ratio depends strongly on the elementary charge 
as ~ and should be significantly suppressed in the case of e* = 2e compared with 

the Fermi-liquid contribution. As a result, the Lorenz number, L = [e/ks)'^ Ke/{Tcr) differs 
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significantly from tlie Sommerfeld value = vr^/S of the standard Fermi-liquid theory, 
if carriers are bipolarons. Here Ke, cr, and e are the electronic thermal conductivity, the 
electrical conductivity, and the elementary charge, respectively. Ref. predicted a rather 
low Lorenz number Lf, for bipolarons, = 6Le/(47r^) ^ O.lSLe, due to the double charge 
on bosonic carriers, and also due to their nearly classical distribution above T^. 

Unfortunately, the extraction of the electron thermal conductivity has proven difficult 
since both the electron term. Kg and the phonon term, Kph are comparable to each other in 
the cuprates. Only recently a new way to determine the Lorenz number has been realized 
by Zhang et al.j2S|, based on the thermal Hall conductivity. The thermal Hall effect allowed 
for an efficient way to separate the phonon heat current even when it is dominant. As a 
result, the "Hall" Lorenz number, Lh = ifi/ksf' tixy/{Ta,j.y)^ has been directly measured 
in Y Ba2Cu-iOQ,<j^ because transverse thermal K^y and electrical axy conductivities involve 
only electrons. Remarkably, the measured value of L^y just above is about the same as 
predicted by the bipolaron model, L^y ~ 0.15Le. 

The experimental L^y showed a strong temperature dependence, which violates the 
Wiedemann-Franz law. This experimental observation is hard to explain in the framework 
of any Fermi-liquid model. Based on the Frohlich- Coulomb model and the conventional 
Boltzmann kinetics in the r approximation we developed a theory of the Lorenz number 
in the cuprates explaining the experimental results by Zhang et al ^29^. We have demon- 
strated that the Wiedemann-Franz law breaks down because of the interference of polaron 
and bipolaron contributions to the heat transport. When thermally excited polarons and 
also triplet pairs are included, our model explains the violation of the Wiedemann-Franz 
law in cuprates and the Hall Lorenz number as seen in the experiment. Fig. 3. 

A large Nernst signal observed well above the resistive critical temperature Tc provides 
another important piece of evidence for a qualitatively different normal state of cuprates 
as compared with conventional superconductors. It has been attributed to a vortex motion 
in a number of cuprates 2^ |^. As a result the magnetic phase diagram of the cuprates 
has been revised with the upper critical field Hc2(T) curve not ending at Tco but at much 
higher temperatures Most surprisingly, Ref. Isil estimated Hc2 at the zero-field resis- 
tive transition temperature of Bi2212, T^o, as high as 50-150 Tesla. However, any scenario 



with a nonzero off-diagonal order parameter (the Bogo^ 



m 



30v-Cor'kov anomalous average). 



F(ri,r2) = ( \E'j(ri)\E'|(r2) ), above Tc, such as of Ref. 
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is difficult to reconcile with 



the extremely sharp resistive and magnetic transitions at Tc in single crystals of cuprates. 
Above Tc the uniform magnetic susceptibility is paramagnetic and the resistivity is perfectly 
normal', showing only a few percent positive or negative magnetoresistivity. Both in-plane 



33, 



35l | and out-of-plane 36| resistive transitions remain sharp in the magnetic field in 



high quality samples providing a reliable determination of a genuine Hc2{T). 

While a significant fraction of research in the field suggests that the superconducting tran- 
sition is only a phase ordering so that a superconducting order parameter F{ri, remains 
nonzero above (resistive) T^, we have recently explained the unusual Nernst signal (which 
is one of the key experiments supporting that viewpoint) as the normal state phenomenon 

El- 

Cuprates are known to be non-stoichiometric compounds. Moreover, undoped cuprates 
are insulators and their superconductivity appears as a result of doping, which inevitably 
introduces additional disorder. Because of these reasons, the theory of doped semiconductors 
might provide an adequate description of the normal state kinetic properties of cuprates. 
Carriers in doped semiconductors and disordered metals occupy states localised by disorder 
and itinerant Bloch-like states. Both types of carriers contribute to the transport properties, 
if the chemical potential fi (or the Fermi level) is close to the energy, where the lowest 
itinerant state appears (i.e. to the mobility edge). Superconducting cuprates are among 
such poor conductors. 

When the chemical potential is near the mobility edge, and the effective mass approx- 
imation is applied, there is no Nernst signal from itinerant carriers alone, because of a 



so-called Sondheimer cancellation 



However, when the localised carriers contribute to 



the longitudinal transport, a finite positive Nernst signal Cy = —Ey/S/^T appears as 

^ = ^r^a,, (5) 
P e 

where p = l/[(2s + l)crsa;] is the resistivity, s is the carrier spin, r is about a constant 
(r^l4.3 for fermions s=l/2, r^2.4 for bosons s=0), and is the Hall angle. Here a^x is 
the conductivity of itinerant carriers above the mobility edge, and ai is the conductivity of 
localised carriers below the edge, which obeys the Mott's law, ai = (Xoexp [— (Tq/T)^/^] (cto 
is about a constant). In two dimensions Tq ~ 8a'^ /{ksNi), where A*"; is the DOS at the Fermi 



level 



3, 



40 



In sufficiently strong magnetic field the radius of a shallow 'impurity' state is about 
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the magnetic length, a ~ (ei?)^/^. Then the normal state Nernst signal is given by 

-^ = a(T)exp[-6(5/T)i/=^], (6) 

where a(T) oc j^^l and h = 2[e/ (ksNi)]^^^ is a constant. As shown in Ref. a 
sing le- p arameter Eq.(6) is in excellent quantitative agreement with the experimental data 
30|,|3j| above the resistive critical temperature Tc{B), Fig. 4. 

As far as statistics of carriers is concerned, the bipolaron or 'preformed boson' picture 
is strongly supported by many other data in underdoped and optimally doped cuprates 
a, while ove.doped cup.ate=^.gM be on the Fe.n,i-h,uM .de. UnUke any fluctuatmg 
preformed pair scenario, eg. |42l|, bosons in our model are perfectly stable, and there is no 
off-diagonal order above their Bose condensation temperature, Tc, There is no off-diagonal 
order parameter above in the overdoped Fermi-liquid either. In both cases the localization 
of carriers by disorder is essential. It is responsible for the strong Nernst signal dependence 
on the magnetic field. Fig. 4. It can be also responsible for a low value of the product 
of the thermopower and the Hall angle S'tanB in some underdoped samples, where the 
contribution to the thermopower from itinerant carriers can be almost cancelled by the 
opposite sign contribution from the localised carriers, so that S oc T at low temperatures. For 
two-dimensional bosons it is feasible even if.„».;0.When.happe.,th:Ne....nal 
is given by Cy = pa^y, where a^y oc r^. Different from that of fermions, the relaxation time 
of bosons is enhanced critically near the Bose- Einstein condensation temperature, Tc{B), 
r oc [T — Tc(i?)]^^/^, as in atomic Bose-gases [i^. Providing StanQn <^ e^, this critical 
enhancement of the relaxation time describes well the temperature dependence of Cy in a 
few Bi2201 close to Tc{B). 

Superconductivity can be readily suppressed by the magnetic field in heavily underdoped 
cuprates to reveal a true normal state down to zero temperature. The normal state in- 
plane resistivity p(T) shows an insulating behaviour at low temperatures, which has been 
understood as the result of elastic (impurity) scattering of nondegenerate carriers above the 
mobility edge [4^. Combining p(T) oc 1/r, Stan 6 oc Tr, and e^-l-S'tanO oc r one obtains 



and 



T 

5'taneoc— — , (7) 



Tq-T 

ey oc (8) 
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where Tq is a fitting parameter. I believe that Eq.(7) and Eq.(8) can account for the exper- 
imental Nernst and 5* tan 6 data in underdoped Lai.92Sro.o8Cu04 and Lai.94Sro.o6Cu04 
I47I 1 with the field- induced insulating behaviour of p{T). The unusual large Nernst sig- 
nal appears due to the absence of the electron-hole symmetry near the mobility edge in the 
transverse transport. Our results strongly support any microscopic theory of cuprates, which 
describes the state above the resistive and magnetic phase transition as perfectly 'normal', 
F(r, r')=0. Unlike 3ll . our model does not require a radical revision of the magnetic phase 
diagram of cuprates j42] . Also an unusual diamagnetic response of underdoped cuprates ex- 
perimentally measured somewhat above Tr(B) 0,0,01 can be understood as the Landau 
diamagnetism of normal state bosons ^511] . 

IV. CHECKERBOARD MODULATIONS OF TUNNELLING DOS 

There are more complicated deviations from the conventional Fermi/BCS-liquid be- 
haviour than the normal state pseudogaps. Recent studies of the gap function revealed 
wo distinctly different gaps with different magnetic field and temperature dependence 



531 154, and the checkerboard spatial modulations of the tunnelling DOS, with |55[ 
and without js^El applied magnetic fields. 

We have proposed a simple phenomenological model |53] explaining two different gaps 
in cuprates. The main assumption, supported by a parameter-free estimate of the Fermi 
energy 5^, is that the attractive potential is large compared with the renormalised Fermi 
energy, so that the ground state is the Bose-Einstein condensate of tightly bound real-space 
pairs. Here I calculate the single particle DOS of strong-coupling (bosonic) superconductors 



by solving the inhomogeneous Bogoliubov-de Genes (BdG) equations [6 

The anomalous Bogoliubov-Gor'kov average F{ri, depends on the relative coordinate 
p = — r2 of two electrons (holes), described by field operators \l/s(r), and on the center- 
of-mass coordinate R = (ri + r2)/2. Its Fourier transform, /(k, K), depends on the relative 
momentum k and on the center-of-mass momentum K. In the BCS theory K =0, and the 
Fourier transform of the order parameter is proportional to the gap in the quasi-particle 
excitation spectrum, /(k, K) oc Ak. Hence the symmetry of the order parameter and the 
symmetry of the gap are the same in the weak-coupling regime. Under the rotation of the 
coordinate system, Ak changes its sign, if the Cooper pairing appears in the d-channel. 
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On the other hand, the symmetry of the order parameter could be different from the 
'internal' symmetry of the pair wave-function, and from the symmetry of a single-particle 
excitation gap in the strong-coupling regime Real-space pairs might have an uncon- 
ventional symmetry due to a specific symmetry of the pairing potential as in the case of 
the Cooper pairs, but in any case the ground state and DOS are homogeneous, if pairs 
are condensed with K = 0. However, if a pair band dispersion has its minima at finite K 
in the center-of-mass Brillouin zone, the Bose condensate is inhomogeneous. In particular, 
the center-of-mass bipolaron energy bands could have their minima at the Brillouin zone 
boundaries at K = (vr, 0) and three other equivalent momenta 0] (here and further I take 
the lattice constant a = 1). These four states are degenerate, so that the condensate wave 
function iplm) in the real (Wannier) space, m = [m^jmy), is their superposition, 

V.(m)= bj^e-^"^"", (9) 

K=(±7r,0),(0,±7r) 

where 6k = ±v^/2 c-numbers, and nc{T) is the atomic density of the Bose-condensate. 
The superposition, Eq.(9), respects the time- reversal and parity symmetries, if 

ilj{m) = ^JtIc [cos(7rm^) ± cos(7rmj,)] . (10) 



The order parameter, Eq.(lO), has (i-wave symmetry changing sign in the real space, when 
the lattice is rotated by 7r/2. This symmetry is entirely due to the pair-band energy disper- 
sion with four minima at K 7^ 0, rather than due a specific pairing potential. It reveals itself 
as a checkerboard modulation of the hole density with two-dimensional patterns, oriented 
along the diagonals. From this insight one can expect a fundamental connection between 
stripes detected by different techniques and the symmetry of the order parameter in 
cuprates |6ll |. 

Now let us take into account that in the superconducting state (T < Tc) single-particle 
excitations interact with the pair condensate via the same attractive potential, which forms 



the pairs 



58l |. The Hamiltonian describing the interaction of excitations with the pair Bose- 



condensate in the Wannier representation is 

H = -J2 - + f^^m,n]clmCsn + [A(m) c{^q^ + H.C.], (11) 

s,m,n m 

where s =t, j is the spin, and A(m) oc ip^m). Applying equations of motion for the 
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Heisenberg operators cJm(^) Csm(^), and the Bogoliubov transformation 6c 

£Tm(t) = 5^[n,(m)a,e-'^''* + vtiin)^^'^'], (12) 

V 



qm(t) = 5^K(m)/5,e-^-* - <(m)ate-^*], (13) 

V 

one obtains BdG equations describing the single-particle excitation spectrum, 

eM(m) = - ^[)f:(m - n) + /i(5ni,n]M(n) + A(m)t;(m), (14) 

n 

-ef (m) = - ^[t(m - n) + /i5m,n]f (n) + A(m)M(m), (15) 

n 

where excitation quantum numbers v are omitted for transparency. Different from the 
conventional BdG equations in the weak-coupling limit, there is virtually no feedback of 
single particle excitations on the off-diagonal potential, A(m), in the strong-coupling regime. 
The number of these excitations is low at temperatures below Ap/A;^, so that the coherent 
potential A(m) is an external (rather than a self-consistent) field, solely determined by the 
pair Bose condensate 



58|. 

While the analytical solution is not possible for any arbitrary off-diagonal interaction 
A(m), one can readily solve the infinite system of discrete equations (14,15) for a periodic 
A(m) with a period commensurate with the lattice constant. For example, 

A(m) = Ac[e^"'"- - e^"""^], (16) 

corresponds to the pair condensate at K = (ivr, 0) and (0, ivr), Eq.(lO), with a temperature 
dependent (coherent) Ac oc ^JrijT). In this case the quasi-momentum k is the proper 
quantum number, i/ = k, and the excitation wave-function is a superposition of plane 

waves, 

«,(m) = Mke''^'" + Mke'^'^-^^-'", (17) 
Here g^, = (7r,0), g^, = (0, tt), and g = (tt, tt) are reciprocal doubled lattice vectors. Substi- 
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tuting Eqs.(17) and (18) into Eqs.(14,15) one obtains four coupled algebraic equations, 

eic^k = CkMk - Ac(vk - Vk), (19) 

Ckitk = Ck-gitk + Ac(vk - Vk), (20) 

-Ckt'k = ^k-g^t'k + Ac(itk - itk), (21) 

-ek^k = Ck-gy^k - Ac(Mk - Mk), (22) 

where ^k = — Z^n^l^)^*"* " ~ A*- '^^^ determinant of the system (19-22) yields the following 
equation for the energy spectrum e: 

(e - Ck)(e - Ck-g)(e + Ck-gJ(e + Ck-gJ 
= A^(2e + Ck-g, + a-gj(2e - Ck - Ck-g)- (23) 

Two positive roots for e describe the single-particle excitation spectrum. Their calculation 
is rather cumbersome, but not in the extreme strong-coupling limit, where the pair binding 
energy 2Ap is large compared with Ac and with the single-particle bandwidth w. The 
chemical potential in this limit is pinned below a single-particle band edge, so that // = 
— {Ap + w/2) is negative, and its magnitude is large compared with A^. Then the right hand 
side in Eq.(23) is a perturbation, and the spectrum is 

A2 

eik ~ Ck - — , (24) 

e2k ~ Ck-g - — , (25) 
A* 

If a metallic tip is placed at the point m above the surface of a sample, the STM current 
I{V, m) creates an electron (or hole) at this point. Applying the Fermi-Dirac golden rule 
and the Bogoliubov transformation. Eq. (12,13), and assuming that the temperature is much 
lower than Ap/ks one readily obtains the tunnelling conductance 

a{V, m) ^ ^^^^ « Yl \M^)MeV - 6.), (26) 

1/ 

which is a local excitation DOS. The solution Eq.(17) leads to a spatially modulated con- 
ductance, 

a{V, m) = areg{V) + amod{V) cos{7Tm^ + TTmy). (27) 
The smooth (regular) contribution is 

areg{V) = (Tq (w^k + «rk)^(eV^ ' ^rk), (28) 



k,r=l,2 
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and the amplitude of the modulated contribution is 

amodiV) = 2ao ^ UrkUrkSieV - erk), (29) 

k,r=l,2 

where ctq is a constant. Conductance modulations reveal a checkerboard pattern, as the 
Bose condensate itself, Eq.(lO), 

a — (Jr 



' reg 



a. 



Acos(7rm^ + irmy), (30) 



reg 



where 



A = 2y^ [uikUiuSieV - eik) + U2uU2kS{eV - esk)] / 



k 



[("ik + Wik)'^(eV^ - eik) + (mL + uli,)6{eV - eak 



is the amplitude of modulations depending on the voltage V and temperature. An analytical 
result can be obtained in the strong-coupling limit with the excitation spectrum given by 
Eqs. (24,25) for the voltage near the threshold, eV ~ Ap. In this case only states near 
bottoms of each excitation band contribute to the integrals in Eq.(30), so that 

^k - eik _ ^ /Qi\ 

Ulk = Uik ^ -Mlk — < ^^Ik, (31) 

eik - 4k-g fJ'W 

and 

^k-g — ^2]i~ ^ ~ Al 

M2k = 1—U2k ~ -M2k < M2k- (32) 

e2k - 4k /iW 

Substituting these expressions into A, Eq.(30), yields in the lowest order of Ac, 

A ^ (33) 

The result, Eq.(30), is reminiscent of STM data where spatial checkerboard 

modulations of a were observed in a few cuprates. Both commensurate and incommensurate 
modulations were found depending on sample composition. In our model the period is 
determined by the center-of mass wave vectors K of the Bose-condensed preformed pairs. 
While the general case has to be solved numerically, the perturbation result, Eq.(30) is 
qualitatively applied for any K at least close to Tc, where the coherent gap is small, if one 
replaces cos(7rma; + 7rmj,) by coslKxirix + Kyiriy) . Different from any other scenario, proposed 
so far, the hole density, which is about twice of the condensate density at low temperatures, is 
spatially modulated with the period determined by the inverse wave vectors corresponding to 
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the center-of-mass pair band-minima. This 'kinetic' interpretation of charge modulations in 
cuprates was originally proposed jfill before STM results became available. It could account 
for those DOS modulations in superconducting samples, which disappear above Tc because 
the coherent gap Ac(T) vanishes, so that A = in Eq.(30). Indeed some inelastic neutron 
scattering experiments show that incommensurate inelastic peaks are observed only in the 
superconducting state of high- Tc cuprates 6^. The vanishing at Tc of incommensurate 
peaks is inconsistent with any other stripe picture, where a characteristic distance needs 
to be observed in the normal state as well. On the other hand some STM studies (see, for 
example 0|) report incommensurate and commensurate DOS modulations somewhat above 
Tc, in particular, in heavily underdoped cuprates [,67]. I believe that those modulations are 
due to a single-particle band structure and impurity states near the top of the valence band 
in doped charge-transfer insulators, rather than a signature of any cooperative phenomenon. 

In conclusion, the strong-coupling Frohlich- Coulomb model, Eq.(l), links charge hetero- 
geneity, pairing, and pseudo-gaps as manifestations of the strong electron-phonon attractive 
interaction in narrow bands of doped Mott-Hubbard insulators. 

The author acknowledges support of this work by the Leverhulme Trust (London) Grant 
No. F/00261/H , The Royal Society, and by EPSRC (UK) Grant No. EP/C518365/1. 
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Figure captions: 



Fig.l. Top panels: small Holstein polaron band dispersions along the main directions 
of the two-dimensional Brillouin zone. Left: A = 1.1, right: A = 1.2. Filled symbols is 
uj = 0.70 1, open symbols is = 0.66 1. Lower panels: the band structure isotope exponent 
for A = 1.1 (left) and A = 1.2 (right). 

Fig. 2. Top panels: small Frohlich polaron band dispersions along the main directions 
of the two-dimensional Brillouin zone. Left: A = 2.5, right: A = 3.0. Filled symbols is 
u = 0.70 1, open symbols is = 0.66 1. Lower panels: the band structure isotope exponent 
for A = 2.5 (left) and A = 3.0 (right). 

n n 

Fig. 3. The Hall Lorenz number Lh |22| fits the experiment (YBa2Cu306.95 |2q|). Charge 
and spin pseudogaps are taken as 675 K and 150 K, respectively, and the ratio of the polaron 
and bipolaron Hall angles is 0.36. The inset gives the ratio of Hall Lorenz number to Lorenz 
number in the model. 

Fig. 4. Eq.(6) fits the experimental signal (symbols) in overdoped Lai.8Sro.2Cu04 [31 1 
with b = 7.32(K/Tesla)^/^. Inset shows a(T) obtained from the fit (dots) together with 
a oc T"^ (line). 
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